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ABSTRACT 

A method which has evolved In our laboratories for the past 20 yr Is 
re-examined with the Intent of Improving Its accuracy and simplicity of 
application to engineering problems. Several modifications are Introduced 
both to the analytical formulation of the Damage Curve Approach, and to the 
procedure for modifying this approach to achieve a Double Linear Damage Rule 
formulation which Immensely simplifies the calculation. Improvements are also 
Introduced In the treatment of mean stress for determining fatigue life of the 
Individual events that enter Into a complex loading history. While the 
procedure Is completely consistent with the results of numerous two-level 
tests that have been conducted on many materials, It Ts still necessary to 
verify applicability to complex loading histories. Caution Is expressed that 
certain phenomena can also Influence the applicability - for example, unusual 
deformation and fracture modes Inherent In complex loading - especially If 
stresses are multlaxlal. Residual stresses at crack tips, and metallurgical 
factors are also Important In creating departures from the cumulative damage 
theories; examples of departures are provided. 

INTRODUCTION 

Treatment of accumulated fatigue damage has received a large amount of 
attention In recent years. The subject has been popular since Palmgren [1] 



first suggested what Is now known as a "Linear Damage Rule." The same rule 
was later Independently proposed by Langer [2] and Miner [3]. However, It was 
soon recognized that while the method has merit for simple treatment of 
complex loading history. Its predictions are often unconservative. Many 
alternative methods of analysis to predict behavior more accurately followed 
and It became necessary periodically to prepare review papers placing all the 
new methods Into perspective, among them the ones by Newmark [4] Kaechele [5], 
Manson [6], O'Neill [7] and Schlve [8]. Schlve's study lists nearly 200 ; 

references pertinent to the treatment of the problem and to the experimental 
programs conducted to evaluate the concepts. 

No comprehensive report has appeared recently to review the considerable 
effort made In the past 12 yr since Schlve's publication. Such a report Is 
long overdue, and It was our original Intent to prepare such a review. 

Because of space limitations, we found, however, that the review had to be 
limited. Thus, In the discussion that follows we shall restrict ourselves to 
elucidation and extension of those aspects of the subject with which we have 
had personal association. Although we will have occasion to refer to the work 
of others, these references will be brief. The more comprehensive review 
remains a subject for the future. 

CUMULATIVE FATIGUE FOR MECHANICAL LOADING 

Most theories of cumulative damage are based on observation of behavior 
at a second loading once the material has undergone loading, short of failure, 
at a previous level. This behavior Is then extended by analogy to treatment 
of a complete history Involving numerous loading levels. In our work we have 
used three basic concepts for the formulation of procedures: the effect of 

prior loading on the convergence and rotation of the fatigue curve In 
subsequent loading, the Damage Curve Approach (DCA), and the Double Linear 
Damage Rule (DLDR) concept. A critique of these methods will now be discussed. 
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Fatigue Curve Convergence and Rotation 

Basis . - Our entry Into the cumulative fatigue damage field arose because 
of observations of a geometric relation between fatigue curves representing 
different degrees of damage either by prior fatigue or by damaging forms of 
machining. Bennett [9] In 1946 published the results of fatigue tests on 
specimens stressed at two successive levels, In Fig. 1. Line PA Is the 
original S-N curve. Specimens prestressed at ±54 k$1 for 33 percent of mean 
life, and then tested at other stress levels, yielded the new S-N curve 
PB. When 90 percent of the life was applied at ±54 ksl, the remaining S-N 
curve was PC. It appeared to us that all three S-N curves projected to a 
single point P at ±60 ksl at 20 000 cycles to failure, and that the concept 
of convergence of S-N curves contained the potential for the formulation of 
a cumulative fatigue damage procedure. 

Another set of curves that suggested to us the convergence concept for 
cumulative fatigue Is shown In Fig. 2 [10]. Here fatigue for a medium 
strength steel are shown for different types of surface finish. As, one 
proceeds from a ground finish through progressively rougher surface finishes, 
the S-N curves simply rotate about a point P. If we regard each form of 
machining as the Imposition of a different level of surface damage similar to 
that which occurs In fatigue, this figure also suggests that each level of 
fatigue damage reflects Its physical effect by causing a rotation of the S-N 
curve for such materials In which the S-N curve can be linearized. 

On the basis of these observations, we proposed an approach to treatment 
of cumulative fatigue damage based on S-N curve rotation [11]. A 
considerable test program was undertaken to check the validity of the 
approach, the major results being shown In Fig. 3. Here, the first loading 
was at ±100 ksl, for which the median life was determined to be 22 000 cycles. 
A series of specimens were therefore subjected to various fractions of this 
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life, and then the S-N curve of the prestressed material was determined for 
each condition of prestress. These S-N curves Indeed appeared to be 
derivable by rotation of the basic S-N curve of the material. Since one 
point on the "remaining" S-N curve Is known by subtraction of the prestress 
cycles from the life at the prestress level, the predicted S-N curve can be 
determined by passing a straight line through this point from the point of 
rotation, P. While this point Is unknown at the outset. It could be 
determined as a material constant, as In Fig. 4, or It could be estimated by 
using the point on the S-N curve at a stress near the ultimate tensile 
strength. Once the rotation point Is known the effect of any loading sequence 
can be determined by successive application of the concept for treating two 
loadings. Each loading uses as Its S-N curve the rotated value already 
produced by the prior loadings, as Illustrated In Fig. 4. 

In Refs. [11] and [12], we assumed that each loading progressively 
reduced the endurance limit. However, for brevity, we shall not review this 
Issue. In many cases the existence of an endurance limit, and Its progressive 
changes. Is not of Importance; however, even If endurance limits need to be 
considered, the rotation concept may not be the best vehicle for doing so. As 
we shall see later, we feel that the convergence framework has other 
limitations as well. Thus we prefer to devote more research to the other 
methods. Interested readers may, however, find Refs. [10] and [11] useful If 
there Is need for Information on changes of endurance limit. 

In a more recent report Hashln and Rotem [13] have presented an 
Interesting discussion of the convergence concept as a special case of a 
framework they have devised for cumulative damage analysis. They avoid the 
Issue of determining the convergence point as a material constant, and take 
this point at a quarter of a cycle. They also speculate on alternative 
formulations, for example, when the convergence point Is at the endurance 
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limit. Fig. 5. Such an approach would not, of course, properly represent the 
usually observed loading order effect (l.e. high stress followed by low stress 
producing a cycle ratio summation less than unity, while a low-high sequence 
produces a sum greater than unity), but It Is an Interesting possibility for 
academic study. Hashln and Rotem also consider the case when stress versus 
cyclic life Is linear on log-log coordinates, rather than semi logarithmic 
coordinates of Figs. 3 and 4. Of course, In this case the constant damage 
curves are presumed straight on the log-log coordinate system. Hashln and 
Rotem do not really require linearity In any particular coordinate system, 
'only that the constant damage curves emanating from the focal point be a 
systematic set, never again Intersecting anywhere. However, If they are not 
linear, then some other criterion must be Introduced to establish the nature 
of the curvature. 

Critique . - Although our Initial efforts at Implementing the convergence 
concept were quite vigorous, several limitations of the concept soon led us to 
seek alternate approaches. One was the obvious contradiction that would occur 
If the prior cycling produced a "remaining life" less than that at the 
convergence point P. For example. In Fig. 6, If a prior loading at the 
stress level of point A were applied to a degree that the remaining life 
were less than 1000 cycles, such as 8, then the damage line after the Initial 
cycling would have a positive slope, BP, which means that higher stress would 
produce longer life - a clear adabsurdum. 

Secondly, as we conducted more tests on other materials It became clear 
that the convergence approach does not always properly represent the damaging 
effect of prior loading on subsequent loading at lower stress levels. We can 
Illustrate this effect by a simple example In Fig. 7. Suppose, for example, 
we apply the Initial loading at the stress level of point A. If we apply 
x percent of the life, at this level, the remaining life, (1 - x) percent at 


5 



this stress level Is shown at point A 1 . If, however, later testing Is at 
stress level B, failure would be expected at B', the linear projection of 
PA 1 to the stress level of B. Experience has Indicated such a test would 
produce a lower life, at B", and In some cases a higher life, as at B"' . 

To Illustrate the somewhat general expectation we resort here to the use 
of the Double Linear Damage Rule which we shall discuss later In the report. 

It Is appropriate to use this rule for Illustration because It reflects more 
accurately the experience we have had with many materials. Figure 8 shows the 
results of detailed computations. Selected life fractions are assumed to be 
applied at stress level A, and the remaining S-N curves are calculated by 
the DLDR approach to be A'B', A"B", etc. Two Interesting results are 
evident: The curves A'B 1 , A"B", etc. are not straight lines, as required by 

the linear-convergence concept; rather they are curved. Also, while A'B 1 
almost converges to P, the curves representing more highly fatigued material 
do not converge. After some fraction of preloading, they start to displace 
somewhat parallel to AB. Hence, the overall effect Is one of rotation plus 
translation. 

Incidentally, the rotation plus translation overcomes the difficulty 
envisioned earlier since the point B can be displaced to any desired life 
level by prior loading without requiring the self-contradictory remaining life 
curve such as BP In Fig. 6. This concept of rotation plus translation was 
first discussed by us In Ref. [14] when we Introduced the DLDR concept. We 
have refined the approach In later publications, as we shall discuss, but the 
basic feature of rotation plus translation Is Inherent In all the variants we 
have studied. 

On the basis of the above limitations of the rotation concept of 
convergence, we have not pursued this .concept further. The alternatives not 
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only give more accurate behavior representation, but are easier to Implement 
as wel 1 . 

Damage Curves 

Basis for single-term damage curves . - The use of damage curves to 
reflect the Interaction of loadings at different stress or strain levels has 
been used since the concepts of cumulative fatigue damages have first been 
considered. In 1948 Rlchart and Newmark [15] Introduced the Idea, but did not 
provide definitive formulae for representing the damage curves In such a way 
as to predict quantitatively the high-low load sequences. Many attempts have 
been made since then to provide suitable analytical relations. In 1981 we 
made our own attempt [16] based on a large amount of experimental data for 
two-level loading conducted In our laboratory. The form we chose, 

D - [da] 

was based on analogy to early crack growth. Figure 9 shows the concept as we 

applied It. The curves represent the damage accumulation for each of the 

2 6 

decade life levels from 10 to 10 . As an example. If cycles are 

3 

applied at the 10 cycles to failure level the damage accumulated will be 

represented by point A. Changing to a second load level, which produces a 

5 

life of 10 cycles to failure, moves us to point B at the same damage as 
A, and cycling can continue for n ? cycles as damage accumulates from B 
to failure at C. Thus, If the damage curve formulae are known, It Is easy to 
calculate n 2 /N 2 once the n^/N^ Is specified If and are known. 

The subscripts refer to sequence of the applied loadings, 1 Is the first and 2 
Is the second level. A more convenient type of plot Is shown In Fig. 10 where 
the life scale Is normalized to failure life. From this figure It Is clear 
that If the first loading Is at the higher stress (lower life) and the second 
Is at the lower stress (higher life) the sum of the cycle ratios will be 


|a Q + (0.18 


V 


IMP - 3 4 5 
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smaller than unity. Since this Is the most commonly observed behavior, the 
qualitative nature of the damage curves Is correct. Before discussing the 
precise quantitative Implications It Is useful to clarify a few points. 

1. Flexibility of constants . - If we express the relation that the damage 
at A and B are equal In Fig. 10, l.e., = D^, then It Is clear that 

the a Q , 0.18 and 2/3 terms In Eq. (1) cancel, and that 



( 2 ) 


More generally for any combination of life levels and we obtain 
the Damage Curve Approach (DCA). 


,0.4 


Hi . i . 

N 2 


(3) 


Thus the relation between n 2 /N 2 and n^/N^ depends only on the ratio 
N.j/N 2 , and Is Independent of any constants that might be Introduced In 
place of the a Q , 0.18, and 2/3. Thus, suppose we replace a Q by zero, 
the 0.18 by 1.0 and the 2/3 by 0.02512 {equal to [1/10 4 ]^' 4 ), the damage 

4 

curve relation for the 10 life becomes linear and the damage curve for any 

other value of N becomes simplified. In fact we can choose to linearize the 

0 4 

damage curve for any by replacing the 2/3 by [1/N^] * . In 

this case 


and 


'ref 


'ref 


(4) 
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Figure 11(a) shows the damage curves for life levels from 10 to 10 7 when the 

4 

damage curve for 10 Is taken as the reference straight line. More 

generally, these same curves are obtained when any life level Is taken 

as the reference life. The damage curves for other lives are labeled as a 

function of N/N^, as seen Figs. 11(a) and (b) which produce the same 

curves for equal values of N/N re f whether the reference life Is taken as 
3 4 

10 or 10 cycles. Only one set of damage curves Is thus needed as long 
as the designation Is according to N/N^. Note that the reference life 
level can be used as the first, second, or any other number In a sequence of 
cumulative loadings. 

2. Significance of the damage curves In relation to two-level loadings . - 
Figure 12(a) shows the damage curves for two-level tests, using the first 
(lower) life ^ as the reference which Is forced to be a straight line. 

If we wanted to use these curves to obtain the n ? /N 2 versus /N^ curves 
for such a test, the procedure would be to go from A to B to produce the 
arbitrary n^N , determine C at the same damage as B, and measure n ? /N 2 
as damage proceeds from C to E. Thus n 2 /N 2 Is equal to length CD. If we 
plot n 2 /N 2 versus n^/N^, as In Fig. 12(b), the point C 1 results. For other 
choices of cycle ratio n^/N^, n 2 /N 2 falls at other points along the curve 
A'C'E'. Examination of Figs. 12(a) and (b) reveals that the curves ACE and 
A'C'E' are similar. One can be obtained from the other by a 90° rotation 
about the central point x. We can then generalize that If the curves of 
Figs. 11(a) and (b) represent the damage line relationships, the expected n -|/ N 1 
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versus n i ^ N 2 curves can obtained by a clockwise 90° rotation. Such a 
family of curves Is shown In Fig. 12(c). 

Basis for a new double-term damage curve system . - Although the 
single-term damage curve system was chosen by us to agree with many 
experimental results generated In two level tests, there Is one feature of the 
system that has been of some concern to us from the very beginning of Its 
use. This feature can be seen In Fig. 12(c). The Initial drop of n 2 /N 2 
for a small value of /N^ Is exstremely rapid. Since the choice of the 
damage curve relation was derived mostly for two-level tests which did not 
Involve such low values of /N^ as to be within the troublesome low 
range, the choice was not In conflict with our experiments. However, after we 
recognized that the resulting feature of the damage curve analysis was 
probably unrealistic, we decided to attempt to Improve Its characteristics. 
Brief attempts to modify the terms In Eq. (1) Indicated that It would be 
difficult to Improve the very low n-|/N-| part of the curve while 
retaining Its good features In the remainder of the range. We decided, 
therefore, that a more appropriate way would be to add a term that would have 
a large significance at low values of n^/N^, but only relatively small 
effect for larger values. In our procedure to alter the Damage Curve equation 
we were guided by our Double-Linear Damage Rule analysis [16] which will be 
discussed later. 

To distinguish the new double-term damage curve from the original 
single-term expression, we refer to It as the Double-Damage Curve Approach 
(DDCA). Details of the derivation are given In Appendix A and only the final 
results are presented here: 


fn_1 

y r 

i Y~1 

fn 1 Y(q 2- 1) 

L N fJ 

+ 

O’ 


l\l 


where 
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Note that q^ Is the slope of the first straight line segment of the D 
versus n/N curve In the DLDR formulation, discussed In the next section. 

The term Is the exponent of the D versus n/N curve In the DCA 
formulation, Eq. (5), and the constant y Is Introduced to represent two 
Intersecting straight lines by a curve having a single methematlcal formula. 

Direct comparison of the newly formulated DDCA, Eq. (6), with the DLDR 

and the original DCA Is given In Fig. 13 for values of N rgf /N of TO -1 , 

-2 -3 

10 , and 10 . As expected, the DDCA equation conforms closely to the 

DLDR In the early portion of the Phase I regime but blends directly Into the 
DCA curves which are also close to the DLDR In Phase II. 

Double Linear Damage Rules 

Background . - One concept that developed early to explain deviations for 
the linear damage rules was that fatigue was at least a two-stage process - 
crack Initiation and crack growth; If the two processes developed at different 
rates for different life levels, there Is no reason for expectation of a 
linear damage accumulation rate for all life levels. An Immediate conjecture 
was that while It Is still possible for linear damage accumulation to develop 
during the Initiation stage, and linear accumulation during the propagation 
stage, propagation might develop at different fractions of life for the 
different life levels. Therefore, a mixture of loadings Involving several 
life levels would not be amenable to analysis by a single linear damage 
analysis. Thus two damage rules were needed; one for Initiation and one for 
propagation. 
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Probably the first to suggest analysis by two linear damage rules was 
Langer [2]. Later, Grover [19] again proposed the concept Independently. 
However, neither proposed quantitative relations for establishing when one 
process ended and the other began, nor how to treat problems Involving more 
than one loading level. Manson [14] made an early attempt In 1966 on the 
basis of limited data, but the relation proposed at that time did not prove to 
be sufficiently accurate to predict results of later experiments, so 
additional attempts [20] were made. Even these attempts were not sufficiently 
general, or of sufficient definitive procedure to be of general engineering 
use. However, In 1981 Manson and Halford [16] refined the formulation and 
procedure to be directly applicable to the analysis of engineering problems. 
Additional valuable contributions both to formulation, and possibly to 
defining the physical mechanisms Involved, were made by Miller and his 
co-workers [21 and 22]. 

Definition of linear and double-linear damage rules . - Rather than 
progressing from damage curve to damage curve as loading Increments 
accumulate, considerable simplification can be Introduced by using the well 
known concepts of single and double-linear damage rules. Consider, for 
example, Fig. 14(a), which Involves two loading life levels. If the damage 
curves are linear, It Is possible to normalize each by dividing by the failure 
life, reducing both damage lines to a single line. Obviously, for this case 
the procedure discussed for summing damage Increments leads to l n/N^. = 

1, l.e., the linear damage rule ( LDR) . Adding loading levels to the history 
does not change this conclusion as long as each damage level Is straight. 

Even If the damage curves consist of broken straight line segments, as shown 
In Fig. 14(b), a linear damage rule will apply If normalizing relative to 
failure life reduces all damage lines to. a single set of broken lines. 

Or, the damage lines are curvilinear over their entire range, as shown 
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In Fig. 14(c), a linear damage rule will apply If all damage lines are reduced 
by normalization to a single curve. Nor does It matter If more loading levels 
are Involved, as long as all reduce to a single curve by one normalization. 

The merit of treating damage accumulation by a linear damage rule rather than 
by damage curves Is great simplification; loading order becomes unimportant, 
and we sum all the cycle ratios without regard to sequence, even when 
Individual loadings at one level are Interspersed with loadings at the other. 
The amount of calculation Involved, and the associated bookkeeping of damage 
accumulation Is, therefore, drastically reduced. 

Now consider the case shown In Fig. 15(a). Two damage curves are shown: 
one Is straight, the other shows a single break. Normalizing by dividing by 
failure life produces the triangular shape In Fig. 15(b). Let us label as 
Dj the region of damage from the origin to the level of AB, and designate 
It as Phase I since the linear nature of OA and OB suggests that In this 
region of damage accumulation the damage curves have something In common. 
Similarly the Phase II region from AB to CE can be labeled Djj. If In 
the Oj region we divide all cycle levels of curve OA by N ft and all 
damage levels by Dj, the line OA" of Fig. 15(c) results. The same 
straight line results when we normalize OB by division by N g and Dj. 

Figure 15(c) now states that as long as the damage has not yet reached 
Dj, we can use a single linear damage rule. In this region. It does not 
matter how many times load level changes from OA to OB, the cycle ratios 
along OA can be summed Independently, and those along OB Independently. 
When the sum of cycle ratios reaches unity, Phase I damage Dj has been 
completed. 

Figure 15(d) Is constructed In a similar manner from the segments AC 
and BE of Fig. 15(a). Cycles are divided by N^ and N gr , 
respectively, and the Increase In damage beyond Dj Is divided by Djj. 
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Since both segments collapse to a single straight line, we can treat damage 
accumulation In this region as independent of application order as well. 

Thus, no matter how many times, or In what sequence, the loading changes from 
the OAC to the OBE level, the damage accumulation depends only on the total 
sum of cycle ratios at each level. But cycle ratios must of course be 
determined relative to the cycles required to apply the damage, not 

total cycles. 

In our early work we designated the Dj region as "crack Initiation" 
and the Djj region as crack propagation. However, we could detect no 
cracks at the transition between the regions, and for other philosophical 
reasons discussed In [16] we have abandoned designations according to physical 
Interpretations, and now simply call these regions I and II. As we shall 
later discuss, the engineering use of the concept Is not restricted by absence 
of physical designation. 

Extension to multiple loading levels . - In the foregoing discussion we 
have limited the consideration to two loads Involving only one break point 
level In the damage relations. As long as breaks occur at the same D level, 
the procedure could be generalized to any number of load levels without change 
In final result. The general case Involving more than two load levels 
requires that the break points be at various levels. Consider first 
Fig. 16(a) In which the damage line OFG has been added. We could, of 
course, divide the damage range Into three regions: below F, between F 

and B, and above B. This would lead to triple-linear damage curves, 

(Fig. 16(b)) and a triple-linear damage rule (Fig. 16(c)). Inclusion of more 
than three loading levels would Involve even a larger number of distinct 
linear regions, and the desired simplicity Is lost. We have, therefore, 
chosen an alternate approach. Only two regions are used, but they are,chosen 
so as to preserve fidelity to the Important damage lines In a given loading 
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sequence and sacrifice fidelity to those loadings which do not have an . 
Important effect. 

The basic philosophy Is Illustrated In Fig. 17. Here we choose OAP as 
the damage line for the life level of the loading which will cause the 

major damage In a particular application. As noted earlier, we can choose 
this damage line as a straight 45° line, ending at point P where D * 1.0 
and n/N 1 = 1.0. The line OBP for a life Is the normalized damage 
line for the loading level that Is of next Importance In producing damage for 
this loading history. (Before performing the analysis the values of N 1 
and may not be known, and an engineering judgment Is made. The 
subsequent analysis provides better choices of N 1 and N^, which quickly 
converge to appropriate values.) As long as only the life levels and 
are Involved In a loading sequence, we can proceed as already 
discussed. If another life level (between and N 2 ) Is added, 
however, Its damage curve will be OCP, with a breakpoint at a higher damage 
level, as previously discussed. One way to keep the breakpoints at the same 
level of damage, Is to replace OCP by OC'P, choosing C 1 where OC 
Intersects the damage level AB. We then have a damage line which coincides 
with the real damage line In region OC' and lies above It In the region 
C'CP. The replacement damage line thus errs on the conservative side since we 
attribute a greater damage to an Increment of loading at this level than would 
occur by using the real damage line. Similarly, If we need to consider a 
loading level with a life greater than we could conservatively replace 
Its damage curve ODP by OD'P. Once we have forced all breakpoints to occur 
at the same damage level we can reduce the problem to one of a double-linear 
damage rule. While the schema of Fig. 17 Is a viable framework for reducing 
any problem to one Involving a double-linear damage rule, the procedure we 
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have actually adopted Is easier to apply and gives very good results, as will 
be discussed. 

Development of Practical Double-Linear Damage Rules 
Earlier work . - Our objective for many years has been to establish a 
practical double linear damage rule framework usable In design and analysis. 
Progress has been gradual. Our first report dates back 20 yr when we first 
proposed [14] that "crack Initiation" and "crack propagation" as the phases 
within which the two linear damage rules applied. Using limited Information 
on one material It first suggested that the crack Initiation period be 

N f - 15N^' 6 whenever > 750 cycles, and equal to 0 for < 750. The 

propagation period would, of course, be 15^" 6 and N^., respectively. This rule 

was simple, but unfortunately not sufficiently accurate for other materials. 

In 1967 a second report [20] examined data on other materials and test 
conditions. We were forced to conclude that not only could we not. In general, 

use the simple formula N p = 15N^' & for the propagation phase, but that It was 

Inappropriate to use the terms "crack Initiation" and "crack propagation" as 
descriptors of the two phases. Not only could we not detect any crack at the 
kneepolnt of a two-level test, when the "crack" should have "Initiated," but 
the kneepolnt shifted as we combined different levels of the second loading 
with Identical first load levels. For example, If a 1000 cycle life level 
loading Is followed by a 10 000 cycle life load level, the kneepolnt Implies 
that "crack Initiation" starts at about 200 cycles. However, If the 1000 
cycle life level Is followed by a loading of life 10^ cycles, the Implied 
"crack Initiation" life for the 1000 cycle level now becomes only about 60 
cycles. Since the material cannot sense while It Is Initially loaded at the 
1000 cycle level what the subsequent loading will be, why should It "crack" 
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after 60 cycles In one case and 200 cycles In the other? For other choices 
of loading levels even greater discrepancies could result. Our conclusion was 
that the type of linearity we perceived for so many materials and loading 
combinations was not a manifestation of "crack Initiation" and "propagation," 
but of some other phenomenon that we did not yet fully understand. 

In the 1967 report we salvaged the concept of a double-linear damage 
rule, but changed the terminology, calling the regions where the two rules 
appeared to be valid as Phase I and Phase II. That we did not know the 
physical significance of these phases did not Invalidate the potential use of 
the concept, but application required an empirical determination of the 
kneepolnt between the two phases for each sequence of loading to be 
evaluated. Unfortunately, such an approach Is not appealing to a designer. 
When one Is trying to choose a material among many candidates, and when the 
loadings may In fact even depend on the material chosen, few designers have 
the Inclination, the facilities, or the lead-time to conduct the required 
tests. 

Recent development . - The matter lay dormant for many years, but Interest 
was again stimulated In this subject during our evaluation of the Damage Curve 
Approach ( DCA) for treating cumulative fatigue damage. The model fashioned 
for DCA was, naturally, also based on the considerable amount of data 
developed for two-level tests on many materials, and the type of results that 
we obtained were as shown in Fig. 18. The continuous curves In this figure 

show the results for the damage curves used at that time. ^We then 

< 

question what these damage curves would look like If we replaced them with two 


+As already discussed, our use of damage curve equations using only a 
single term gave somewhat unrealistic results for low values of n-j /N-j . 

Use of the double damage curve equation (6) discussed earlier produces results 
very similar to the dotted lines In this figure. 
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straight lines, and sketched the dotted lines shown In the figure. Of special 
Interest was that the damage curves for a sequence Involving two load levels 
N.j and depended only on the ratio N^/^, and not on the 
Individual values of and • The same curve would be obtained If 
the first life level were 1000 cycles and the second 10 000 cycles as would 
result when a 100 000 cycle loading was followed by a 10 6 cycle loading, 
therefore we would have to expect the same Implication for the DLDR. 

This observation was precisely the clue we needed to understand the 

confusing results derived by stating a DLDR based on crack Initiation and 

crack propagation concepts. Referring to the Illustration used earlier. If a 

1000 life level loading were followed by a 10 000 cycle level, the N^/N^ 

ratio Is 0.1, and we should expect the apparent knee to be at about 200 

cycles. But If the 1000 cycle life loading level Is followed by a 10 6 life 

-3 

level, the N^/N 2 ratio becomes 10 , and we should expect the apparent 

knee to be at about 60 cycles. It was not a matter of "Initiating" a crack at 
either of the two knee points. Rather, It was a question of which N-j/f^ 
curve was Involved. By defining Phase I as the region up to the kneepolnt, 
and Phase II as the region beyond the kneepolnt, we can avoid the 
philosophical question of Its physical significance and yet devise a simple 
but useful analytical procedure. 

We returned to our raw data to see If an equation could be derived for 
the coordinates of the kneepolnt based on the parameter N -|/ N 2 ra ^er 
than on the N 1 parameter as we had previously sought. The results are 
shown In Fig. 10. Indeed the data of the kneepolnt coordinates normalized 
better bn the basis of the ratio N^/N^. (By comparison we show In 
Fig. 20 that when only was used to correlate the horizontal coordinate 
of the kneepolnt, as required If crack Initiation occurs here, the results are 
poor.) 
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From Fig. 19 the coordinates of the kneepolnt are 



Once the kneepolnt Is established, It becomes possible to re-examine all the 
two-level tests we have conducted to determine If the data could have been 
predicted In advance. Figure 21 shows how double-linear relations, using 
Eq. (7) to establish the kneepolnt, have worked for three materials - 300CVM, 
SAE4130, and T1-6A1-4V. The agreement generally Is quite good. 

Application to multiple loadings . - We next consider how to treat 
multiple loading levels. 

1. Multiple changes In loading Involving two load levels : Consider a 
case Involving a mixture of loadings which If applied Individually would 
produce 1000 and 100 000 cycles to failure. In this mixture changes occur 
from one load to the other many times In arbitrary order. A simple case, for 
example, Is one In which blocks consisting of 1 percent each of the basic 
loadings are successively applied until failure occurs. Each block contains 
10 cycles of the loading which alone produces a 1000 cycle life, and 1000 
cycles of the loading which alone produces 100 000 cycle life. We seek to 
determine how many blocks can be sustained. Figure 22 shows how this problem 
would be organized. 

The first step Is to reinterpret the DLDR In terms of equivalent damage 
curves. Figure 22(a) shows the double-linear representation on the basis that 

3 

Phase I Is consumed entirely at the 10 cycle-to-fal lure level, and Phase II 

5 

at the 10 cycle-to-fallure level. The coordinates of point B are 

computed by Eq. (7), and the numerical values are shown In the lower part of 

the figure. Figure 22(b) shows the equivalent damage lines obtained by a 90° 

counterclockwise rotation of Fig. 22(a), as previously discussed. Thus 

3 

A'B"C' represents damage accumulation at the 10 life loading, while 
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A'B'C shows how damage accumulates for the 10 life loading for a problem 
which contains only these two loadings. 

The next step Is to track the damage accumulation as successive loading 

blocks are applied. In the early loading the damage accumulates along A'B" 

3 5 

when the 10 loading Is applied, and along A'B' when the 10 loading Is 

applied. Thus the straightforward way Is to move between these two lines as 

block after block Is applied, until the damage reaches the level B''B'. 

Beyond this point the pertinent damage lines to use are B"C' and B'C*. 

Again the straightforward procedure Is to move back and forth between these 

two lines. Increasing damage until point C Is reached. Proceeding In this 

way requires many changes from one line to the next (over 52 times). 

However, If we recognize that until the damage reaches B " B ' , both damage 

lines are straight, then a linear damage rule Is valid and the loading order 

does not matter. We can then assume that the same effect will be produced by 

3 

separately lumping together all the cycles of the 10 life level and all the 
cycles of the 10 5 level. Thus If It requires B^ blocks to complete 

3 

Phase I, then 10 B^ cycles of the 10 life level and 1000 cycles 

5 

of the 10 life level will be applied. Since the number of cycles to 

complete Phase I for the 10 3 life level Is 111, and for the 10 5 life 

level. Is 79 445, the equation for B^ In Fig. 22 applies, and It Is seen 
+ 

that 9.74 blocks are required to complete Phase I. By similar reasoning, 

It follows from the calculations of Fig. 22(b) that 16.67 blocks are required 
to complete Phase II, for a total number of 26.41 blocks. This Is In contrast 
to 50 blocks that would be unconservatlvely predicted by the Miner Linear 
Damage Rule (LOR). 


+For simplicity we accept decimal values of blocks, even though a 
nonintegral value means that some loadings of one level are required to 
complete the block without being matched by the other loading. 
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The simplicity of this type of approach Is evident for the elementary 
problem Illustrated here. When extended to more complex loading histories the 
saving of calculation labor Is even more striking. 

2. Loadings Involving three load levels : In Fig. 23 we have added a 
third load level which by Itself would produce an Intermediate life of 10 000 
cycles, and again assume block loadings Involving 1 percent of Individual 
lives. As seen In the figure, B 1 and B" are obtained as before. The 

4 

problem Is with the damage curve that Is to represent the 10 cycle life 

3 

level. Its kneepolnt, based on Individual combination with the 10 cycle 
life loading Is at D, which does not have the same vertical coordinate as 
B'B". To use a DLDR we must find a suitable replacement for the line ADC 
having a break-point along B'B". Several alternatives are possible. One Is 
to choose D 1 at the average where AD and CD Intersect the B'B" line so 
that AD'C replaces ADC. A second way Is to use the Interpolation formula 
shown In Fig. 23 which was derived In [16]. Basically, It Is an analytical 
relation for the cycle ratio to the end of Phase I which Is consistent with 
the formulas used to get the coordinates of B' and B" and which both 
Interpolates and extrapolates In a manner consistent with experimental 
experience. This matter will be discussed In connection with Fig. 24. 

For the case of the 10 000 cycle life In Fig. 23, the formula produces a 
location for D 1 at n/N^ = 0.476 which Is nearly the same as averaging 
D" and D‘". The replacement damage line AD'C In Fig. 23 Is seen to be 
below the real damage line ADC In the Phase I range, and above It In the 
Phase II range, compensating for errors Introduced by forcing all breakpoints 
to have the same ordinate. 

Once we have forced all the breakpoints to be at the same ordinate, 
linear damage rules can be expressed for the Phase I region below this 
ordinate, and for the Phase II region above It. Damage calculations are shown 
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In Fig. 23 for the case of blocks of 10, 100 and 1000 cycles applied 

3 4 5 

respectively at 10 , 10 , and 10 life levels. As before, we first find 
out how much Phase I damage a single block of loading produces, and then 
determine how many blocks It takes to complete Phase I. Similarly we 
calculate how many blocks are required to complete Phase II, and therefore the 
total number of blocks to failure. In this case, as seen In Fig. 23, a total 
of 20.7 blocks are required. The LOR would have predicted 33.3 blocks, about 
60 percent higher than the DL0R computation. To treat this simple problem by 
the DCA would have required changing progressively from one to the other of 
the three damage curves, approximately 60 times. Even for this simple case 
the calculations would Involve an order of magnitude more labor than the DLDR 
procedure. But the result would have been nearly Identical to that obtained 
by the DLDR. For the DCA analysis the prediction Is 21.0 blocks; for the DDCA 
the prediction Is 23.3 blocks, and for the DLDR analysis It Is 20.7 blocks. 

3. Cases Involving large number of loading levels : In the more complex 

problem Involving many loading levels two Important questions arise: 1) which 

two load levels should be used as the baseline values around which to 
establish when Phase I ends and Phase II begins, and 2) how to construct the 
damage curves for loading levels In the extrapolated range, that Is lower than 
the lower life of the reference damage lines or higher than the higher life of 
the reference damage lines. 

Construction of damage curves . - We address the second of these questions 
first by extending the simple problem we have been discussing using four 
loading levels shown In Fig. 24, loadings which alone would produce lives of 
10 3 , 10 4 , 10^, and 10^ cycles. Also, to avoid the Issue of which 
loading produces most of the damage, we assume that each block contains 
1 percent of the failure cycles for each loading level Individually. 

3 

Figure 24(a) shows how each of the damage lines would be drawn If the 10 
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damage line Is taken as the 45° reference line. Note that the kneepolnts do 
not have the same value, as follows from Eq. (7). One of these kneepolnts 
needs to be chosen, thus allowing us to retain only two damage curves, but the 
other two must be compromised. 

In Fig. 24(b) the kneepolnt associated with the 10 3 and 10 & life 

4 5 

levels Is retained, and the 10 and 10 life damage lines are compromised. 

They are constructed In the same manner as used In Fig. 23, using the formula 

shown In this figure. The calculations for blocks to complete Phases I and II 

are then made In an Identical manner as In Fig. 23, resulting In 11.47 blocks. 

In Fig. 24(c) the 10 3 life damage line Is paired with the 10 5 life 

4 6 

damage curve, thus requiring the compromising of the 10 and 10 damage 
curves, - one In the Interpolated range of the faithful damage curves, and one 
In the extrapolated range. Again the same formulas shown In Fig. 23 are used 
to establish the location of the compromised damage curves, and the procedure 
for determining number of blocks Is the same. For this case the computed 
number of blocks Is 12.03. 

Finally, In Fig. 24(d), the 10 3 and 10 4 life damage lines are 

5 6 

faithfully retained, and damage curves for 10 and 10 life are 
compromised. Again, of course, the formulas of Fig. 23 are used to establish 
the coordinates of the kneepolnts, and the computation procedure Is the same 
as In the previous cases. Half of the damage curves are now In the 
extrapolated range, but the computed number of 13.77 blocks Is not far 
different from the other. 

It should be recognized that this example represents an extreme case. 

Here the contribution of each loading to the damage accumulation Is the same, 
since In each block all contribute equal numerical values of damage. Yet, 
even though equally damaging damage lines were compromised, the computed 
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blocks differed by less than 10 percent. But the conclusion Is clear that If 
the reference damage curves are at least as Important In contributing damage 
as those damage lines that have been compromised, the use of our basic 
equation shown In Fig. 23 Is adequate for establishing the damage lines used 
In the analysis. 

Choice of the reference damage lines . - It stands to reason that the 
degree of error Introduced by compromising fidelity of a damage line depends 
on how Importantly the loading associated with that damage line contributes to 
the total damage. Thus the most Important loading events should be used as 
the reference uncompromised damage lines while the lesser events can be 
compromised without Introducing much error. If, as seen In the example of 
Fig. 24, compromising Important lines do not seriously affect the resulting 
calculations, then compromising the lesser contributors should not be too 
serious. But first we have to know which are the Important events. It would 
seem that the answer could be obtained by a successive approximation 
approach. The first step could be an ordinary linear damage analysis, and the 
cycle ratio at failure for each of the contribution events determined. The 
two most damaging events can be taken as and N 2> where Is the 
lower of the two. If there Is ambiguity because several events contribute 
approximately the same damage, then choose and as far apart as 
possible, so that other Important damaging events will have life values 
between them. Then make a DLDR analysis, and review which events are now the 
most damaging. If other loadings now appear to be better candidates as N.j 
and N 2> redo the analysis. Since each analysis Is quite simple, consisting 
of only two sets of computations each of which Is the same as a LOR analysis, 
this procedure should not be too demanding. Further refinement could be added 
by making successive Iterations until life predictions change little and the 
final choice of and Involves the most damaging events or unless 
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they are recognized as the proper choices according to the engineering 
Judgment of the analyst. 

4. Complex service cycle : In Ref. [16] we analyzed a cycle that had 

originally been discussed In Ref. [23] representing a complex mission cycle 
for a jet engine. The data for the calculations are shown In Table I. Every 
block of loading Involved 14 events, each of known total stralnrange Ac and 
tensile mean stress Cg, columns (2) and (3). Each event was applied a 
known number of times, per mission column (4). From these values, life was 
calculated* for each event, column (5). The procedure used then was to 
choose the cyclic life of the event that produced the lowest life as the 
N-| value, In this case event No. 8 with a life of 2500 cycles. Similarly, 
the value was the highest one for the events, 64 000 cycles for event 
6. From these two values, N ^ and Nj ^ for each event was calculated 
from the equations shown In Fig. 23. The remaining procedure was Identical to 
that described In earlier paragraphs. The n/N^ and n/N^ values for 
each event were calculated, columns (8) and (9), leading to a calculation of 
the number of missions required to complete Phase I and Phase II, 
respectively. Thus, for this problem 79 missions were required to complete 
Phase I and 200 missions for Phase II, leading to a total of 279 missions to 
produce failure. 

It Is now appropriate to carry the calculations one step further, and 
revise the procedure to determine and N 2 more appropriately. A 
better choice of and N 2 can be made by considering the higher and 


+In our earlier analysis we used a different formula for life when mean 
stress Is present. Now that we have revised the procedure to Include mean 
stress effects, as discussed In Appendix B, the values In column (5) would be 
somewhat different. However, since the changes would not seriously affect the 
numerical values Involves, and would have no significance on procedure, we 
have not remade the calculations for column (5). 
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lower of the two life values for the events that produce the largest amounts 
of damage, rather than the events that produce the highest and lowest cyclic 
lives . 

Using the rationale already discussed, we have added In Table I the 
column (10) whlch’ls not contained In Ref. [23]. After the first calculation 
Is made, as earlier proposed, and the number of missions required to complete 
each of the two phases are calculated, the damage for each type of event Is 
calculated (column (10)). Events 4 and 8 are the most damaging. Hence we 
remade the calculations for = 2500 and = 5550. On this basis, 

events 4 and 8 are still found to be the two most damaging, thus Iteration of 
calculation Is not necessary. The total number of missions now calculated Is 
277, which differs little from the 279 missions earlier determined, and gives 
us confidence that we can err considerably on the choice of optimum values of 
N-j and N^ without affecting the final life calculation significantly. 
However, the degree of error may depend significantly on the type of complex 
cycle analyzed. For example. If the complex cycling Involves a preponderance 
of damage accumulation at a single life level, errors due to choice of 
and N 2 , or even of choice of damage rule, will become Insignificant. 

Impetus to study cumulative damage rules for complex loading Is provided 
by the extreme high- and low-cycle fatigue operating conditions experienced by 
components In the NASA Space Shuttle Main Engines (SSME). Several thermal 
low-cycle fatigue and mechanically Induced high-frequency, high-cycle fatigue 
are superimposed on turbine blades and other critical structural components. 
One such component, called a L0X Post, Is a slender tube carrying cryogenic 
oxygen through Its bore while exposed to high-temperature, hydrogen-rich steam 
on Its exterior. Flow-Induced vibration provides the high-cycle fatigue 
excitation while the thermal transients of engine firing and shut-down 
contribute to the low-cycle fatigue loadings. To understand better the damage 
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Interaction In the alloys of construction, cumulative damage experiments were 
conducted on both 316 SS and a cobalt-base superalloy, Haynes 188. Results 
are presented In Refs. [25] to [27]. In view of the Improvements (DDCA) to 
the OCA, we have re-examined both sets of cumulative damage data and have 
sought optimum material specific constants. Figure 25 shows the results for 
316 SS at 1300 °F compared with calculations by the DDCA, Eq. (6). The 
exponents a and 8 from Eq. (6) have been optimized for the data shown, and 
have values of 0.23 and 0.63, respectively. Agreement Is better than obtained 
earlier using the original DCA formulation and the single constant, 8 = 0.4. 
The added flexibility of another constant In the DDCA promotes Improved 
accuracy of representation. 

Similar treatment, Ref. [27], was given to the analysis of data on Haynes 
188 reported by Blzon, et al. [26]. The optimum constants are a = 0.35 and 
8 = 0.60 In contrast to the universalized values of 0.25 and 0.40, 
respectively. Comparison of the DDCA calculations with the experimental 
results are presented In Fig. 26. A band of expected behavior Is shown, 
reflecting the fact that tests were performed at somewhat different values of 
N^/N . Again, good agreement Is obtained. Most Importantly Is the 
unconservative deviation (up to a factor of 10 In cyclic life) of the 
experimental results from the classical LDR. The nonlinear features of the 
DDCA (and the DIOR) accurately model the Interactions of high- and low-cycle 
fatigue. 

DISCUSSION 

In this re-examination of our procedures for computing damage 
accumulation we have Introduced three major changes. 

Restating the DCA by adding a term which enables us to make the damage 
curve nearly Identical with the DLDR that has been found characteristic for a 
large number of materials In two-level fatigue tests. The new equation Is 
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called the Double Damage Curve Approach (DDCA). We have also defined more 
clearly the relation between damage curves and their approximation by the DLDR 
Introducing a successive-approximation procedure for the choice of the 
parameters of Phases I and II of the DLDR. 

Modifying the procedure for Including mean stress effects In cyclic 
life. A few comments are appropriate for each of these Items. 

Restatement of Damage Curve Approach and Relation to DLDR 
It Is Important to emphasize that we treat each event on the basis of 
life level, not the parameters that lead to Its life level. The strain range 
and maximum, minimum, or mean stress are Included In the determination of the 
life level of the event. Thus any appropriate theory can be applied to 
combined them to determine their life level. Cumulative damage Is then 
calculated In the same manner as It would be at the same life level produced 
by other parameters leading to this life. This procedure Is In contrast to 
other approaches, e.g. Dubuc, et al. [28] and Bui Quoc [29, 30], which 
Incorporate the complex stress and strain parameters Into the expression for 
the damage function. The advantage of our procedure Is that the approach may 
change for determining the life level without affecting the damage curve 
analysis. For example, as already discussed, we have progressed to a more 
accurate and generalized method of treating mean stress, but the change In 
Itself does not alter the DCA. 

However, we changed the functional form of the damage curve 
representation to make It conform more closely to extensive data, and to make 
It more consistent with the DLDR. The new functional form Is called the 
Double Damage Curve Approach (DDCA). Thus when we make a comparison between 
results obtained by the DDCA and DLDR analyses we do not have to sort out how 
much of the difference Is due to the simplifications used to establish the 


28 


DLDR, and how much /Is due to the difference Implicit In the two types of 
damage function used. 

The change In our new function Is consistent with some of the thinking 
that has recently been Introduced by other Investigators. For example. Miller 
and Zacharlah [21] and Miller and Ibrahim [22] devised damage functions which 
are also based on two stages of damage accumulation - development of shallow 
surface cracks and propagation of a macroscopic dominant crack leading to 
failure. In their first effort, damage curves were represented by straight 
lines on log-log coordinates, with breakpoints arranged In similar fashion to 
ours, as seen In Figs. 23 and 24; that Is, the line joining the breakpoints 
had a negative slope. However, In their later work they changed their damage 
functions so that the slope of the line joining the breakpoints Is positive. 

It Is Interesting that while the results of damage curve analyses are similar 
when based on either our curves or those of Miller and his co-workers, the 
governing parameters are different. If we use crack length as the measure of 
damage our crack lengths become macroscopic fairly early, while Miller and 
co-workers derive crack lengths of ultra-microscopic size In the Phase I 
stage, and macroscopic size only very close to failure. The Interesting point 
Is that, as discussed earlier In the report, the damage accumulation process 
In going from one damaging curve to another does not depend on all the 
parameters defining the damage curves, but only on one or two critical 
constants. Thus they can derive the same final result event though the 
baseline Information differs appreciably from ours. 

Successive Approximation Procedures for Choice of and N 2 

In our previous publications we recognized the Importance of considering 
the lowest and highest lives to establish the values of and In 

the DLDR formulation. A restriction was Invoked requiring 1 percent life 
fraction be Imposed before a level could be considered an extreme. We have 
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since come to realize that life levels for which the greatest amounts of 
damage (life fraction, n/N) are accumulated are also of Importance because 
they represent the major participants In the problem. Thus, we have suggested 
an Iterative process that homes In on the two most damaging life levels and 
uses them as the and life levels In the DLDR. Calculations 
presented In this paper have borne out the accuracy of this approach. 

Modifying Procedure for Including Mean Stress Effects 
The cumulative fatigue damage approaches proposed In this paper do not 
rely upon how mean stresses affect cyclic life, or the specific mean stress 
formulation employed. Nonetheless, a cumulative damage analysis cannot 

proceed until all the life levels can be quantified. Since mean stresses 

\ 

strongly Influence fatigue life, a reliable, yet general, mean stress 
formulation Is a necessary adjunct to the cumulative damage approaches. We 
have modified the mean stress approach from what was used In our earlier 
publications, and have arrived at a more general means for representing these 
Important effects. It has been demonstrated In Appendix B that a wide range 
of classical mean stress effects can be accommodated by the newly adopted 
approach. Because of the approach's flexibility, however, It does require a 
small amount of fatigue data Involving mean stress to evaluate the constants. 
Evaluation techniques are discussed In Appendix B. 

CAUTIONS 

While we have been active In developing a procedure for treating damage 
accumulations either through damage curves or through a DLDR based on these 
curves, few complex service histories have been treated by these methods. 

Much experience should be gained In diverse applications before these or other 
such methods, can be accepted for general use. In many cases, It can be 
expected that the method will be successful by the very fortuitous nature of 
such complex loading histories. They can contain many events, but In fact 
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most will contribute negligibly to the damage, and life will be governed by 
one or two of the major loading events. The results of the Teledyne study (In 
Table I) are typical of what can happen In such cases. Current methods are 
satisfactory for such cases since they were developed from two load level 
tests. 

The real caution comes In connection with circumstances that have not yet 
been factored Into this or other methods. A partial list of such 
possibilities Is: 

Effect of Stress Multlaxlal 1 t y 

This subject Is currently receiving Intensive study that Is beyond the 
scope of this paper, and may reveal unusual and d 1 f f Icult-to-predlct results, 
especially when loadings In one set of directions or slip systems are 
subsequently followed by loadings In other directions favoring other slip and 
fracture systems. 

Deformation Systems In Complex Loading Not Revealed In Individual Loadings 

Wood and Relmann [31] tested copper In torsion wherein large loading 
amplitudes for a portion of the life were then followed by lower loading 
amplitudes. The result did not follow a simpler linear damage rule, 

Fig. 27(a). Relmann [32] later conducted similar torsion tests of Iron. More 
cycles at the lower strain level were possible after the higher strain level 
had been applied than If It had not been applied. Fig. 27(b). We have tested 
copper and steels under similar conditions In axial loadings, but did not 
reproduce the type of effect observed In torsion. Wood explained his results 
on the basis that when the larger amplitudes were first applied the strain 
tended to be absorbed by many slip-planes. Once activated to slip, they 
remained effective In slipping when the lower amplitude was later applied, 
even though these planes might not have participated If the lower amplitude 
loading were applied singly. Thus the plastic strain of the lower loading was 
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distributed over more slip planes when preceded by the higher loading, 
decreasing the strain per slip plane, and Increasing the life. Why we have 
not seen such a phenomenon In the many materials we have tested axially Is not 
clear. 

Residual Stresses 

The Importance of residual stress has been especially revealed through 
the recent emphasis on the fracture mechanics of cracked structures. One such 
experience Is demonstrated by Fig. 28. Two-level load tests were conducted on 
aluminum In bending [33]. We Introduced a notch on only one side of the 
specimen, and ran the high load first to the point at which a small crack had 
developed at the root of the notch. The last of the high level loadings 
applied was with the notch In tension, producing a residual compressive stress 
at the tip of the crack when the load was removed. When the lower load was 
applied, the remaining cycle ratio was not less than unity, as normally 
expected, but It was greater by more than 100 than the Initial life of the 
material at the lower load. The crack had been arrested by the residual 
compression, and did not progress when the lower load level was applied. The 
failure at the low load Initiated at the opposite side of the specimen, when 
there was no notch, and progressed toward the crack that had been started by 
the high load. The DLDR certainly was not directly applicable here. 

Metallurgical Instabilities 

Stress and strain together with exposure to high temperature reactive 
gases can produce metallurgical Instabilities such as strain aging. Some 
materials, subjected to strain at high temperature, develop precipitates on 
the dislocations that alter the subsequent deformation and strength 
characteristics of the material. Thus loadings which Induce such 
precipitations during one portion of the loading history will reveal a 
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behavior In another portion which by Itself would be different because of no 
metallurgical precipitation. An example Is presented below. 

While Investigation creep-fatigue Interaction In the cobalt-base 
superalloy L-605, at 1400 °F,.we observed a strong strain-aging effect [34]. 
Carbon, Initially In solution with the matrix, precipitated forming an 
M„C, type carbide along the generated dislocations. As the carbides 

do b 

Increased In number and size during cycling, the cyclic stress-strain response 
of the alloy changed dramatically. The alloy was thus Increasing Its cyclic 
flow resistance and altering Its resistance to fatigue failure. The amount of 
change depended upon several factors - amount and duration of deformation and 
temperature. The Implication to cumulative damage analysis Is that the 
fatigue life at a prescribed load or strain level Is not a unique quantity 
dictated by the magnitude of the loading. Instead, details of prior loadings 
can alter the fatigue life relation. For example, consider two high 
temperature loading levels. One under high strains gives rise to dislocation 
generation and carbide precipitation, producing a life, N 1 . The other Is at 
such a low strain level that no precipitation nuclei are formed, the material 
does not harden and the ensuing life Is N > . If a few high strain 

cycles are Imposed on another sample, carbides will form, the material 
hardens, and Its fatigue resistance will be altered. When the partially 
fatigued sample Is subsequently loaded at the smaller strain level. Its 
fatigue resistance should no longer be associated with N ? . Competing 
processes are occurring simultaneously; fatigue crack nucleatlon and growth 
(damage) and material hardening. The difficulty to the analyst Is In how to 
separate these factors to accurately predict remaining fatigue life. Further 
research remains as to how best to resolve such complex problems. 

Another example of altered fatigue resistance at high temperatures Is one 
encountered while studying Stralnrange Partitioning for creep-fatigue 
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analysis, Ref. [35]. When 316 SS was tested with hold periods at peak 
compressive stress, an oxide surface layer formed, which cracked during the 
tensile portion of each cycle. After a few cycles of loading, the surface was 
riddled with short cracks, most of which hadn't as yet penetrated the 
substrate. For small cyclic strains, oxides still form, but would not crack 
during the tensile excursions. Hence, a sample loaded at large strains for a 
fraction of Its life and then fatigued to failure at a small strain range 
would experience an additional damaging phenomenon (a cracked oxide surface 
layer) not accounted for In the original fatigue curve of the alloy at the low 
strain level. Determining the correct life level for cumulative damage 
analysis could become a tedious task. Our approach does not specifically 
address the Issue raised above. Instead, the approach Is more general and 
deals only with life levels and not with how the life level Is attained. 

The phasing of temperature and strain cycling during thermal fatigue can 
produce additional complications to the problems of cumulative damage 
assessment. In a study, Ref. [36], of thermomechanical fatigue behavior of 
the nickel-base superalloy, MAR M 200, we observed significant differences In 
life depending upon whether In- or out-of-phase cycling was Imposed. During 
In-phase cycling, early developed cracks became filled with oxide while being 
held open by the tensile stress. This situation Is not experienced In the 
out-of-phase cycling for which compressive stresses close cracks at the high 
temperature, retarding oxidation within the cracks. Out-of-phase cycling was 
considerably less damaging than In-phase. Furthermore, the two phasings 
result In different modes of cracking at the ubiquitous Internal carbides. 
In-phase cycling produced carbide-matrix Interfaclal cracking while 
out-of-phase cycling caused the carbides themselves to fracture. The 
Implication of the above findings to cumulative damage analysis are of 
concern. For example, consider two straining levels, Ae^ and Ae^, 
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giving rise to the same cyclic life, N f . One condition Involves In-phase 
cycling, and the other, out-of-phase, Ae 2 [Ac 2 > Ae^]. If cycling of a 

sample Is started at A c 1 and Is discontinued after a life fraction /N^ and 

Is resumed at Ae 2 until failure occurs, the remaining life fraction according 
to the current rules (and the LDR) would simply be n 2 /N^ = 1 - n^/Np l.e., 
n 1 + n 2 = N f Since both conditions produce the same life, any combination 
would be predicted also to produce the same life. However, It Is unlikely 
that such a simple result would be borne out by experiment, due to different 
damage accumulation mechanisms existing for the two loading conditions. With 
different mechanisms of damage, the damage curves of D versus cycle fraction 
n/Nf would not coincide for the two conditions, and hence damage summations 
would differ from 1.0. 

The lesson to be learned from these examples Is that the cumulative 
damage rules developed In this paper which are based upon a single basic 
fatigue crack Initiation and propagation mechanism, will require refinement 
for application to unique high temperature conditions wherein additional 
damage mechanisms can come Into play. 

CONCLUDING REMARKS 

In this report we have described a philosophy of analysis that has 
evolved In our laboratories over more than 20 yr. Our emphasis has been 
simplicity of engineering application, and minimizing the baseline Information 
required for Implementation. While much of the methodology has been discussed 
In other reports, this presentation takes the evolution process a few steps 
further. The format of the damage curve analysis has been altered and the 
procedure clarified to be consistent with the Implied damage curves used In 
the Double-Linear Damage Rule (DLDR). We call the new formulation the Double 
Damage Curve Approach (DDCA). We have gained Insight as to how to Iterate the 
choice of these substitute damage curves to minimize error while still 
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retaining simplicity In analysis. Since our basic damage curves are now 
consistent In both the ODCA and OLDR analyses, the usefulness of the damage 
curve analysis procedure has been Increased. By making the same calculations 
through the OOCA we can establish how much error Is Introduced by the 
compromises brought about to alter them for DLDR analysis. On the other hand, 
It also becomes clear from such calculations how much simplicity Is gained by 
the OLDR procedure compared to retention of damage curves with attendant 
computational complexity. 

An Important feature of our methods Is the characterization of an event 

only by Its life, not the parameters that enter Into determining the life. 

Thus a stralnrange and associated mean stress that leads to a given life Is 
treated In Identical manner as a smaller stralnrange with associated larger 
tensile mean stress which also leads to the same life. In this way 
Improvements In life calculation per se do not later the equations Involved In 
the DOCA or DLDR analyses. We have. In fact, discussed what we regard to be 
Improvement In accounting for mean stress. Thus, while the life values that 
enter Into a damage calculation may be altered by using the new mean stress 
relations, the equations operating on these life values do not change. 

Experience with the application of the method to a spectrum of complex 

loading types Is currently limited, but we hope to make detailed computations 

In generic cases to evaluate the effects of various parameters. Sample 
computations are Included In this report. 

Finally, we urge caution In the use of the method described, or Indeed 
any other method, to Insure that some unexpected phenomenon, not Inherently 
contained In the framework of the method, Is Introduced Inadvertently. Among 
these are stress multlaxlal 1 ty , unidentified deformation and fracture 
mechanisms, unknown residual stresses (especially at notches and crack tips), 
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a Mean stress. 




APPENDIX A 


DERIVATION OF THE DOUBLE DAMAGE CURVE APPROACH ( DDCA) 

Our choice of terms to enter the two term damage curve approach was such 
as to produce a continuous curve coincident with the DLDR at small values of 
life fraction and coincident with the single term DCA equation at large life 
fractions. The first term Is linear In n/N and the second Is a power 
function of n/N 



(A-l ) 


where: 



The coefficient Is the slope of the first damage accumulation line In 
the DLDR. The second term In Eq. (A-l) Is the same as the original single 
damage curve term except for a reduction of Its coefficient from 1 to 
(1 - Q-j) to force conformity with D = 1 at n/N = 1. N ref Is the 
reference life condition for which damage could be considered to be 
accumulating linearly toward unity, and N Is the life level of Interest. 

While Eq. (A-l) serves to Improve the slope of the damage curve at Its 
origin, especially for large values of N/N re ^ (Fig. A-l(a)), additional 
Improvement Is required for lower values of N/N ref (Fig. A-l(b)). In both 
cases, the double term equation deviates the greatest amount In the transition 
region between the two terms. 
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To achieve Improved fit over the entire range we resorted to a previously 
used approach [17] and [18]. Given an equation where the dependent variable 
Is expressed as the sum of two power laws of the Independent variable, 

y = X m + X n ( A-6) 

It Is clear that y -* X m when X n « X m , and y -» X n when X m « X n . Thus y 
becomes nearly coincident with X m at one extreme of X and nearly coincident 
with X n at the other extreme. Plotted on log-log coordinates y becomes 
nearly straight In two regions of X. Between the two extreme regions y 
changes slope gradually from m to n. However, we can force y to conform 
to X m and X n over a wider range, and change more abruptly between the two 
slopes by rewriting the equation In the form 

y T = X 1 ™ * X yn <A-7> 

when X m » X n , y still approaches X m , and similarly when X n » X m , y ■* X n . 
The asymptotic lines are thus preserved. However, by the proper choice of y 
the curve of y versus X can be altered In the range where X m s x n . 

Figure A-2 shows the application of this concept for the case at hand. 

Rewriting Eq. (A-5) , 


D 


n 


q{ ♦ [i 




T[q 2 -1] 




( A- 8 ) 


If we choose several values of y and plot D versus n/N for N/N^ = 10, 
we get the curves shown In Fig. A-2. It Is clear that y = 5 gives a 
sufficiently close fit to the double linear damage line, and we have, there- 
fore, tentatively settled on this value. 
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APPENDIX B 


MODIFIED PROCEDURE FOR MEAN STRESS EFFECTS 
Service cycles Involving variable amplitude loading will likely encounter 
mean stresses during portions of the loadings. Since mean stresses 
significantly alter the expected fatigue life relative to a zero mean stress 
condition, techniques must be Implemented to determine the associated life 
levels for use with the cumulative damage approaches described In the body of 
this report. The following describes a promising modified procedure for more 
reliably dealing with mean stress effects on fatigue life. Consider a 
material for which the Manson-Cof f In-Basquln diagram has already been 
established, as shown In Fig. B-l(a). No mean stresses are Involved In this 
figure. Figure B-l(b) Is exactly the same plot, except that the life scale 
has changed to reflect the presence of mean stress In accordance with the 
formulation discussed by 

A+BlogN 

Heldmann [24]. The multiplier on the life scale Is [1 - (o Q /o^) r ]. 

Although o ^/ E and b are known from Fig. B-l(a) the constants A and B 
are at the outset unknown. Several tests Involving mean stresses are 
necessary to determine these constants. As an example, assume these tests are 
conducted with a mean stress ratio of 0.2, l.e., o /c, = 0.2, o = 0.2 x 

o T o 

130 = 26 ks 1 . For Illustration assume that two tests are conducted, one with 
an alternating stress of amplitude 65 ksl, the other with an alternating 
stress amplitude of 40 ksl. The lives are measured, and tabulated as shown In 
Table B-l . We treat three cases In which the measured lives are assumed to be 
those In Table B-l. Consider the analysis of the datum for Case I where 
La / 2 = 1 30( 2N^ ) “^ . 1 0 _ ^ ^ y^e elastic strain amplitude for the 

3 

alternating stress Is 65/30x10 = 0.002167 In. /In. We thus proceed to 

point P In Fig. B-l(b), Indicating that the numerical value on the 
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horizontal axis Is 1024. Thus we apply this value and the measured life 
(2Nj. = 2 x 55) to the horizontal scale 

1024 = (2 x 55) [1 - (0.2) A+B1 ° 9(55) ] ^ 1Q) (B-1) 

from which we calculate 

A + 1.74B =1.00 (B-2) 

If we apply the same procedure to the second datum for Case I, where 
Ao/2 = 40 ksl, we develop the equation 

A + ( 3 . 85) B = 1.00 ( B-3) 

Solving for A and B from Eqs. (B-2) and (B-3). 

A = 1.0; B = 0 

In a similar manner. If we treat the two measured points for Case II, we get 

A = 2.0; B = 0 

and finally for Case III, the measured life values give 

A = 3.0; B = -0.42 

Case I, therefore, represents a material of the Goodman-Morrow type. All the 

normalized alternating stress/mean stress diagrams consist of a single 

straight line, as In Fig. B-2. Case II represents a material for which the 

normalized diagram consists of a single curve of the Gerber type, with convex 

curvature. Finally, Case III represents a material of general behavior, 

wherein the normalized alternating stress/mean stress diagram consists of a 

family of curves, one for each life level. Once the material has been 

characterized, any corresponding type of plot can be made that suits the needs 

of the Investigator. Figure B-3 shows the mean stress effect for each decade 

3 6 

life level between 10 to 10 cycles. 

Figure B-1 can also be used to analyze data wherein the stralnrange and 
mean stress are known. Thus If point Q Is known from a total stralnrange 
(which may Include appreciable plasticity), we can still establish the 
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abscissa at point R. If mean stress Is also known, together with the values 
of A and B for an already characterized material, N^. can be determined 
from a simple transcendental equation. Alternatively, If life Is known from 
characterization tests for which Ac and a Q are also known, values of A 
and B can be computed In the same manner as already Illustrated. 

Obviously, for the schematic Illustrations, only two tests are needed for 
determining A and B. In practice, many tests may be used to get A and 
B from a least squares solution. Or, In fact, It may develop that these many 
tests will define a more general function for N^. Instead of the relation 
A + Blog(N f ).. 
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TABLE B-l. - HYPOTHETICAL FATIGUE LIVES UNDER 
VARIOUS MEANS STRESS CONDITIONS. EXAMPLES 
USED TO EVALUATE CONSTANTS A . AND B IN 
GENERALIZED MEAN STRESS EQUATION 



Nfo 

Nf, 

cycle life 

°a. 

ksi 

o 

II 

o 

o 

o 0 = 26 ksi 

( 

Case I 

Case II Case III 


+65 

512 



+40 

65 736 



Computed 



A 


1.00 

2.00 3.00 

B 


0 

0 -.42 














, ALTERNATING STRESS AMPLITUDE. 




Fig. 1. Fatigue curves for specimens pre- 
fatigued at ± 54 ksi for 0, 33, and 903b of 
crack initiation life. After Bennett (9). 


Fig. 2. Surface finish effect on fatigue of 
medium strength steel. After Noll and 
Erickson (10). 



Fig. 3. Comparison of S-log N lines for original and pre-stressed SAE 
4130 heat treated material. Cyclic pre-stress conditions were 100 ksi. 
Median points shown in all cases. After Manson, Nachtigall and 
Freche(ll). 



, ALTERNATING STRESS AMPLITUDE 



Fig. 4. General location of S-log N lines for 
original material and materials subjected to a 
single and double pre-stress condition after 
Manson, Nachtigall, and Freche (11). 



N, CYCLES TO FAILURE (LOG SCALE) 


Fig. 5. Damage lines converging at the endurance limit. After Hashin 
and Rotem(13). 


P (N = 1000 CYCLES) 



N, CYCLES TO FAILURE 

Fig. 6. Example of physically unacceptable damage line obtained by 
rotation about convergent point only. 


, ALTERNATING STRESS AMPLITUDE, ksi 



Fig. 7. Convergence approach over-predicts remaining life at 
stress level B. 



Fig. 8. Rotation and displacement of fatigue damage lines calculated using the double linear damage 
rule. 








(b) One combination of Nj and N 2 . 

Fig. 12. Damage curves and damage interaction. 




n 2 /N 2 , REMAINING LIFE FRACTION 



n^Nj, INITIAL LIFE FRACTION 
(c) All combinations of Nj and N 2 . 
Fig. 12. Concluded. 



Fig. 13. Double damage curve approach (DDCA) blends the 
DLDR at low cycle ratios with the DCA at high cycle ratios. 












, DAMAGE 



N a N c N b N e 0 N A /N C Ng/Ng 1 ' 0 


N, CYCLES CYCLE FRACTION 

(a) Damage accumulation curves. (b> Double-linear damage rule lead- 

ing to non-linear life fraction 
summations. 

(END OF PHASE I) (END OF PHASE II) 

A", B" C", E" 



PHASE I CYCLE FRACTION PHASE II CYCLE FRACTION 

(c) Linear life fraction rule for (d) Linear life fraction rule for 

phase I damage accumula- phase II damage accumula 

tion. tion. 

Fig. 15. Elements of double-linear damage accumulation. 



(a) Damage accumulation curves. 


Fig. 16. Elements of triple-linear damage accumulation. 






PHASE I DAMAGE FRACTION 


C'.E'.G' 



(b) Triple-linear damage curves leading to 
non-linear life fraction summations. 


Fig. 16. Continued. 
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(c) Triple-linear damage rule. 
Fig. 16. Concluded. 
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Fig. 17. Conservative approximation leading to a double-linear damage rule 
when breakpoints in damage curves are not at the same level. 
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APPLIED CYCLE RATIO. nj/Nj 

Fig. 18. Cycle ratio relationship for two level 
tests as deduced from damage curve approach, 
and replacement of damage curves by two 
straight line segments to create a double linear 
damage rule. 
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(a) Exponent. (b) Coefficient. 

Fig. 19. Determination of exponent and coefficient in the equations used to calculate the kneepoint for the 
double linear damage rule. 



LIFE AT FIRST LOADING LEVEL, Nj 

Fig. 20. Lack of correlation between kneepoint 
coordinates and the concept of a discrete crack 
initiation event. 
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(a) Maraging 300 CVM steel. 

Fig. 21. Two load level test results for three engineering alloys showing comparison with predictions by the damage curve approach and 
the double linear damage rule. Data from Manson, Freche, and Ensign (20) and previously unpublished NASA results. 
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SAE 4130 STEEL (SOFT) 
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A ROOM TEMPERATURE 

O 77 °C (170 °F) 

DAMAGE CURVE ANALYSIS 

DOUBLE LINEAR DAMAGE RULE 



(b) SAE 4130 steel. 


(c) Titanium alloy 6AI-4V. ^ 


Fig. 21. Concluded. 



(a) Conventional interaction lines. (b) Reinterpretation as damage lines. 

N i - lOOO; » n - 111; Hj, l - 881 

N = 100,000; Nj » 79,445; ? » 20,555 

Applying 1 % of single loading life in blocks to failure: 

Bj (10/111) ♦ (1000/79,445) “ 1.0; B r = 9.74 Blocks 

B (10/B89) ** (1000/20,555) = 1.0; B n » 16.67 Blocks 

Based on DLDR, Total Blocks at Failure = Bj + Bjj " 26.41 
Based on Miner Linear Damage Rule, Total Blocks at Failure = 50. 

Fig. 22. Double linear damage rule applied to block loading involving two loading levels. 
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(10/111) +(100/4,900) + (1000/79,445) = 1.0; Bj = 8.13 Blocks 

B n (iO/8 89) +(100/5,092) + (1000/20,555) * 1.0; Bjj «*12.57 Blocks 

Based on DLDR, Total Blocks at Failure = B^ + B^ = 20.7 
Based on Miner Linear Damage Rule, Total Blocks at Failure «* 33.3 

Fig. 23. Double linear damage rule applied to block loading involving three loading levels. 
Third level intermediate to nominal extreme life levels Nj and 1^. 
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n/N f , CYCLE FRACTION 

(a) Representation of double linear damage rule by damage lines. 

Fig. 24. Double linear damage rule applied to block loading involving four loading levels. 
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(b) Application of double linear 
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TOTAL LIFE 
N f 

PHASE I LIFE 
N I 

PHASE II LIFE 
N II 

10 3 

62 

938 

10 4 

3745 

6255 

10 5 

70658 

29342 

10 6 

884412 

115588 


Calculation of Number of Blocks if Apply 196 
Life Fraction at Each Life Level Shown: 

For Phase I: 

(10/62) + (10 2 /3745) + (10 3 /70658) + (10 4 /884412) = 

0.2135 

Thus, 4,68 Blocks for Phase I. 
for Phase II: 3 4 

(10/938) + (10 2 /62S5) + (10 /29342) ♦ (10 /115588) = 

0.1472 

Thus , 6.79 Blocks for Phase II 
TOTAL BLOCKS = 11.47 


damage rule with Nj ■ 10^ and N 2 ■ 10^. 


TOTAL LIFE 
N f 

PHASE I LIFE 
N i 

PHASE II LIFE 
N II 

3 



10 

111 

889 

10 4 

4910 

5090 

5 



10 

79445 

20555 

6 



10 

928300 

71700 


Calculation of Number of Blocks if Apply 1% 
Life Fraction at Each Life Level Shown: 


For Phase I : _ 3 u 

(10/111) + (10^A910) + (10 /7944S) + (10 /928300) = 

0.1338 

Thus, 7.47 Blocks for Phase I 
For Phase II : _ u 

(lo/to) + (io z /5090) + (10V20S55) + (10V71700) = 

0.2191 

Thus, 4.56 Blocks for Phase II 
TOTAL BLOCKS = 12.03 


(c) Application of double linear damage rule with Nj • 10^ and ■ 10 5 . 



TOTAL LIFE 
N f 

PHASE I LIFE 
N ! 

PHASE II LIFE 
N 

II 

10 3 

197 

803 

4 

10 

6345 

3655 

10 5 

88044 

11956 

6 

10 

964987 

35013 


Calculation of Number of Blocks if Apply 196 
Life Fraction at Each Life Level Shown: 


(10/197) + (10 2 /6345) + (10 3 /B8044) + (10 4 /964987) « 

0.0882 

Thus, 11.33 Blocks for Phase I 
For Phase II :•» 3 4 

(10/803) + (10 /3655) + (10 /119S6) + (10 /35013) = 

0.4091 

Thus, 2.44 Blocks for Phase II 
TOTAL BLOCKS = 13.77 


(d) Application of double linear damage rule with Nj - 10 3 and N 2 * 10 4 . 
Fig. 24. Concluded. 


n 2 /N 2 (HCF) n 2 /N 2 , SECOND CYCLE FRACTION (HCF) 



Fig. 25. Improved representation of HCF/LCF inter- 
actions using DDCA with optimized constants 
(“ - 0.23, 0 * 0.63, y = 5.0). Comparison with LDR 
and DCA (0 = 0.4). 



Fig. 26. Correlation of HCF/LCF interaction by the 
double curve approach (DDCA). Data from ref. (26). 
Cross-hatched area represents range of N}/N 2 
data. 
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(a) Results for copper, after Wood and 
Reimann (31). 
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lb) Results for iron, after Reimann (32). 

Fig. 27. Unusual cumulative fatigue damage 
behavior in torsion. 



(a) 65-20 ksi. (b) 40-20. (c) 65-40. 


Fig. 28. Residual stress effects with 7075-T6 aluminum, ref. (33). 









Fig. A-2. ODCA, Equation (A-2) follows DLDR for low 
values of damage and follows DCA for high values of 
damage for y • 5. 



(a) Zero mean stress condition. 



(b) General mean stress condition. 

Fig. 8-1. Manson-Coffin-Basquin fatigue relation. 




Fig. B-2. Normalized mean stress behavior for cases 
I and II. Valid for all life levels. 



o 0 /o { 

Fig. B-3. Normalized mean stress behavior for case 
III. Shape of curve changes with fatigue life 
level. 
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